Chaos in effective classical and quantum dynamics 
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We investigate the dynamics of classical and quantum TV-component 0'' oscillators in the presence 
of an external field. In the large A'^ limit the effective dynamics is described by two-degree-of-freedom 
classical Hamiltonian systems. In the classical model we observe chaotic orbits for any value of the 
external field, while in the quantum case chaos is strongly suppressed. A simple explanation of this 
behaviour is found in the change in the structure of the orbits induced by quantum corrections. 
Consistently with Heisenberg's principle, quantum fluctuations are forced away from zero, removing 
in the effective quantum dynamics a hyperbolic fixed point that is a major source of chaos in the 
classical model. 
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0^ The study of the quantum mechanics of those systems whose classical counterpart exhibits chaotic dynamics has 

attracted a lot of interest in recent years, and is an open and rapidly evolving field |l]j|]. Chaos does not exist in the 
linear evolution of the quantum state vector, hence different approaches to identify quantum features that correspond 
to classical chaos have been developed, ranging from the application of random matrix theory to the statistical analysis 
of energy spectra ||^ and to various semiclassical approximations Q . Dynamical chaos in the actual quantum evolution 
' may show up in mean-field approaches [^, or using Bohm's formulation of quantum mechanics to define quantum 
trajectories and quantum Lyapunov exponents ||^. Moreover, at a semiclassical level, the dynamics of quantum 
expectation values can be chaotic |@-^. 

In this Letter we consider the effective dynamics of quantum expectation values as obtained in the large N limit. 
' The purpose of the present work is twofold. First, we want to compare the quantum effective dynamics of a model 
system with the classical effective dynamics of the same system at the same level of approximation, in order to 
C ■ unambigously detect the effect of the quantum corrections on dynamical chaos. This effect turns out to be a strong 
suppression of chaos with respect to the classical case. Second, we want to show that in our model the suppression of 
chaos in the quantum dynamics has a clear physical origin in the fact that quantum fluctuations must be nonvanishing, 
i.e., in the Heisenberg principle, 
'"pi ■ As a model system we consider a A^-component field theory in d -|- 1 space-time dimensions, in the presence of 
' , an external field B. We shall limit ourselves to the case d = 0, which allows us to offer a simple intuitive explanation 
of the effectiveness of quantum corrections in suppressing chaos. The Lagrangian that we consider is (a = 1, . . . , iV) 

> 

X ... 

We perform a 1/A^ expansion, keeping only the leading order term in both the classical and quantum case. In the 
. - - former case we start by writing 

= ^0/3 + = + (50a (2) 

and we approximate the quadratic fluctuations by considering all of them equivalent in the large N limit 

Hi^e- (3) 

By inserting Eqs. (||) and (||) in Eq. (|l|) it turns out that the dynamics of the mean field Lp and of its root mean square 
fluctuation ^ is governed by the following effective Hamiltonian 

n^\{.' + V^) + ^{^'+e-vlf-B^. (4) 



where the two canonically conjugated pairs of variables are tp, tt and ry, and vq = jX is the minimum of the 

potential energy in Eq. for _B = 0. 

In the quantum case we consider the time evolution of the expectation value = '^fy.{4>a) /N from a given initial 
quantum state. This initial value quantum problem can be formulated by using the "closed time path" functional 
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formalism . The application of this formalism to the present case was developed by Cooper et al. , who showed 
that the evolution equations in the large N limit are (classical) Hamilton's equations for the effective Hamiltonian 
(we keep the dependence on the external source B), 

n=\{.^+v') + l{v^+e-vlf + ^-B^^ (5) 

where ct = 2n + 1, n = (a^a) being the expectation value, in the initial state, of the particle number operator for a 
single oscillator, and ^ the expectation value of the r.m.s. fluctuation of the fields, in close analogy to the classical 
case. In Ref. it was also shown that Ti in Eq. (^ is just the expectation value of the full quantum Hamiltonian 
in a general mixed (initial) state characterized by a Gaussian density matrix. 

We notice that at this order of the expansion the correspondence between the classical and the quantum case is very 
strict, since both are described in terms of corresponding degrees of freedom, the mean field Lp and its r.m.s. fluctuations 
^, and the classical effective Hamiltonian is just the quantum one with S = 0. This shows that the approximation (^ is 
equivalent to retaining only the classical contribution to the fluctuations. The quantum correction to the Hamiltonian 
(^ keeps the fluctuation ^ away from zero, consistently with Heisenberg's uncertainty principle. In the following we 
will always refer to the Hamiltonian (||), distinguishing between the classical and the quantum case according to the 
value of h. 

In the case of a vanishing external field {B = 0) the Hamiltonian (^ is integrable in both the classical {h = 0) and 
quantum case [h^Q). The integrals of motion are the total energy E and the function 

I={^^-i^f + ^^^ . (6) 

The conservation of / in Eq. (|^) is due to the fact that the Hamiltonian (|^) can be seen as the Hamiltonian, 
in cylindrical coordinates, of a particle in three dimensions moving in a central potential, being the azimuthal 
coordinate and ^ the radial one. The term Y^a^ jAS^pP' is indeed the centrifugal barrier term which shows up 
passing from cartesian to cylindrical variables. From conservation of angular momentum it follows that 1 = — L\ 
is conserved. 

Switching on the external field (5^0), the two systems become nonintegrable. They have two integrable limits. 
These limits are: harmonic oscillators as £' — » 0; and the integrable Hamiltonian with i? = as £^ ^ cxo. In fact, 
as the perturbation is linear in ip^ it will become negligible with respect to the other terms in Ji for sufficiently large 
E. Hence, we expect that chaotic orbits may show up in an intermediate energy range, whose width will depend 
on B. We have studied the dynamics by numerically integrating the canonical equations of motion derived from the 
Hamiltonian (|^) using a bilateral symplectic algorithm [ p2| . The values of the parameters were fixed by working in 
natural units h = \ (in the quantum case) and by chosing \ — vq = a = \. The choice a = 1 corresponds to an initial 
vacuum state for the number operator a^a []lT| . 

We have studied chaos from both a qualitative and a quantitative point of view, i.e., we have calculated Poincare 
sections jisj in the plane (^, 77) and measured the Lyapunov exponent of every single orbit. This has allowed us to 
also obtain an estimate of the relative measure fi of the chaotic regions in the phase space, defined as the ratio of the 
number of trajectories whose Lyapunov exponent is positive to the total number of trajectories, so that < < 1. 
For each value of the energy E and of the field B we have estimated n{E, B) from a sample of 1000 orbits picked up 
at random in the allowed region on the section surface. We have considered the energy range Q < E < 1, always fixing 
the zero of £' = as the minimum energy allowed in the classical i? = case, and field intensity range < -B < 0.5. 

In Fig. |l| we plot ^{E, B) for the classical and the quantum case. It is evident that, as soon as i? 7^ 0, chaotic orbits 
suddenly appear in the classical case. For small values of i?, i.e., B — 0.01, such orbits are present only in a small 
interval of energies centered around E ~ 0.1, then the chaotic energy interval broadens as B grows, and eventually 
fills the whole explored energy range as B becomes larger than 0.3. In the quantum case the situation is completely 
different: no chaotic orbits are detected for B < 0.3, then chaos appears but at considerably larger values of E as 
compared to the classical case. To give an example, in Figures I and I we show a comparison between Poincare 
sections in the classical and the quantum case at two different values of the external field (notice that the external 
field B does not affect the shape of the region of the plane (C,??) accessible to the system, since B is coupled to (p). 

In order to quantify the degree of chaos at a given energy E and external field B we have considered also the 
ensemble average (A)(£', B) of the Lyapunov exponent over the samples of 1000 trajectories used to compute fi{E, B). 
A comparison between the classical and the quantum case is reported in Fig. ^ for the same values of B as in Figs. 
^ and ^. Looking at these figure it is evident that chaos is strongly suppressed in the quantum case with respect to 
the classical case. 
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As i? ^ oo the quantum and classical models are equivalent, hence it is worth considering, in addition to the 
average measures of chaos B) and {X){E, B), also an average measure of the relative importance of the quantum 
part of the Hamiltonian, given by the ensemble average {Q){E, B), where for each single orbit Q is defined as 

Here Vq — {h'^a'^)/ {8£,'^) is the quantum correction to the potential, and V is the total potential — suitably normalized 
| p5[ in order that < Q < 1 — and (•)( stands for a time average along the orbit. The parameter Q{E, B) has a 
smooth dependence on both E and B, at variance with the other parameters (Fig. ||). The transition from completely 
ordered, to mixed (ordered + chaotic), to almost completely chaotic dynamics that is observed in the quantum case 
as B, E > 0.3 does not correspond to any transition from mainly quantum to mainly classical dynamics. 

We now give a simple and intuitive explanation for the suppression of chaos by the quantum correction, in our 
model. Let us consider the map of the plane (^, 77) obtained by a Poincare section of the Hamiltonian system defined by 
Eq. (H), such as those reported in Figs. || and ||. As i? = both the classical and the quantum system are integrable. 
Hence the trajectories of the map lie on invariant tori. Nevertheless, the geometry of such tori is dramatically different 
in the two cases: in the classical one, as E > 1/8 a hyperbolic fixed point at X = (0,0) exists, and the trajectory 
that passes through X is actually a separatrix. Such a hyperbolic point is due to the presence of a local maximum in 
the potential. The quantum correction to the Hamiltonian is a "centrifugal term" that removes the local maximum 
of the potential and replaces it with an infinite barrier. Consequently, no hyperbolic fixed point in (0,0) exists in the 
quantum Poincare section. As soon as the perturbation —Bip is turned on, chaos immediately shows up in the classical 
case just in the neighborhood of X, because the stable and unstable manifolds, which constituted the separatrix of the 
unperturbed map, split and have infinite intersections [l^ ]. In the quantum case this major source of chaos is removed 
because no separatrix exists in the unperturbed case, and chaos shows up only when the perturbation has completely 
distorted the original shape of the potential. In physical terms, the quantum suppression of classical chaos in our 
model is due to the fact that the quantum fluctuations are kept away from zero — consistently with Heisenberg's 
uncertainty principle — by a quantum term in the effective Hamiltonian. 

To summarize, we have presented an example in which the phenomenon of quantum smoothing of classical chaos 
not only clearly shows up, but also finds a simple explanation on physical grounds. In order to understand how 
general this explanation can be, further work is needed. On one hand, the large N expansion belongs to semiclassical 
approximations, the validity of which, as far as chaos is concerned, can be questionable. In fact, in some cases it 
has been explicitly found that the onset of chaos is in correspondence to the breakdown of the approximation, and 
that the exact evolution of the quantum expectation values is not sensitive to the initial conditions |^,^. Yet it 
has been argued that "semiquantum chaos" can be a real effect in open quantum systems since they are driven in a 
semiclassical regime by the interactions with the environment P,p7| . On the other hand, our results on quantum and 
classical Lyapunov exponents (see e.g. Fig. ^ are in qualitative agreement with those reported in Ref. where a 
quantum Lyapunov exponent — defined via the Bohm approach to quantum mechanics — was found positive but 
smaller than the classical one in a model of an hydrogen atom in an oscillating electric field. In Bohm's theory jlj] 
particles obey classical equations of motion with an additional force derived from a "quantum potential" that is of 
order as the quantum correction to the effective Hamiltonian (^ is. Bohm's equations of motion are exact, being 
completely equivalent to the standard quantum theory, but to write these equations for the system (|l|) would require 
the solution of the full time-dependent Schrodinger equation, thus the analysis of the exact Bohmian dynamics of our 
model system is practically unfeasible. Our treatment is approximate but tractable, and the features it shares with 
Bohmian mechanics are certainly suggestive. 
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FIG. 1. Relative measure /j, of the chaotic component of phase space vs energy E at different values of the external field B. 
Each point is an average over a sample of 1000 randomly chosen orbits, (a) classical case; (6) quantum case. Symbols in both 
cases: B = 0.01 (circles), B — 0.05 (solid circles), B = 0.1 (triangles), B = 0.3 (solid triangles), B = 0.5 (squares). Errorbars 
are of the same size as the data points. 
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FIG. 2. Poincare section of the Hamiltonian flow defined by Eq. (g) witfi B = 0.05. (a) classical case, (b) quantum case; the 
energy is i5 = 0.25 in both cases. 



0.5 1 1.5 




0.5 1 1.5 3 



FIG. 3. The same as in Fig. | with B = 0.5 and E = 0.6. 
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FIG. 4. Comparison between the values of {X){E) in the classical (open circles) and quantum (solid circles) case at (a) 
B = 0.05 and (6) B = 0.5. 
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FIG. 5. Quantum parameter (Q) (see text) vs the energy E at different values of the field B. Symbols as in Fig. Inset: 
(Q) vs the field B for two values oi E, E = 0.2 (solid circles) and E = 0.8 (circles). 



6 



